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Abstract 

The effect of subcritical liadron bubbles on a first-order quark-hadron phase 
transition is studied. These subcritical hadron bubbles are created due to ther- 
mal fluctuations, and can introduce a finite amount of phase mixing (quark 
phase mixed with hadron phase) even at and above the critical tempera- 
ture. For reasonable choices of surface tension and correlation length, as ob- 
tained from the lattice QCD calculations, we show that the amount of phase 
mixing at the critical temperature remains below the percolation threshold. 
Thus, as the system cools below the critical temperature, the transition pro- 
ceeds through the nucleation of critical-size hadron bubbles from a metastable 
quark-gluon phase (QGP), within an inhomogeneous background populated 
by an equilibrium distribution of subcritical hadron bubbles. The inhomo- 
geneity of the medium results in a substantial reduction of the nucleation 
barrier for critical bubbles. Using the corrected nucleation barrier, we es- 
timate the amount of supercooling for different parameters controlling the 
phase transition, and briefly discuss its implications to cosmology and heavy- 
ion collisions. 

PACS number(s): 12.38.Mh, 64.60.Qb, 05.70.Fh, 25.75-q, 98.80.Cq 
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I. INTRODUCTION 



The hadronization of Quark Gluon Plasma (QGP) possibly produced in the early uni- 
verse or expected to be formed in relativistic heavy-ion collisions has been the focus of much 
attention during the past few years. However, the mechanism of hadronization (QCD phase 
transition) remains an open question. The prediction of lattice QCD on the order of the 
transition is still unclear, if physical masses for quarks are used Quenched QCD (no 
dynamical quarks) shows a first-order phase transition, albeit a weak one, with small surface 
tension and latent heat 0. Assuming the transition to be first-order, homogeneous nucle- 
ation theory |^,^ has been invoked extensively to study the dynamics of the quark hadron 
phase transition both in the context of early universe as well as for the plasma produced 
during relativistic heavy- ion collisions 0-0. In this picture, the transition is initiated by 
the nucleation of critical-size hadron bubbles from a supercooled metastable QGP phase. 
These hadron bubbles can grow against surface tension, converting the QGP phase into the 
hadron phase as the temperature drops below the critical temperature, Tc- This is indeed 
the case for a sufficiently strong first order transition, where the assumption of a homoge- 
neous background of QGP is justified at the time when the nucleation begins. However, 
for a weak enough transition, the QGP phase may not remain in a pure homogeneous state 
even at T = Tc, due to pre-transitional phenomena. For temperatures much above Tc, 
matter is in a pure QGP phase with the effective potential exhibiting one minimum at 0=0. 
Here is an effective scalar order parameter generally used to model the effective potential 
describing the dynamics of a phase transition. As the plasma expands and cools to some 
temperature Ti, an inflection point is developed away from the origin which on further cool- 
ing separates into a maximum at 4>=(j)m and a local minimum at = 0/i, corresponding 
to the hadron phase. At T = Tc, the potential is degenerate with a barrier separating 
the two phases. "Pre-transitional phenomena" refers to nonperturbative dynamical effects 
above Tc in the range Tc < T < Ti. Such phenomena are known to occur in several ar- 
eas of condensed matter physics, as in the case of isotropic to nematic phase transition in 
liquid crystals [1^, and are also expected in the cosmological electroweak phase transition 



leading to large phase mixing at T = T^; WM. In such cases, the phase transition may 



proceed either through percolation [T^IT^ or, if the phase mixing is below the percolation 
threshold, by the nucleation of critical bubbles in the background of isolated hadronic do- 
mains, which grow as T drops below Tc- In either case, the kinetics is quite different from 
what is expected on the basis of homogeneous nucleation |T^. We will argue that, for a 
wide range of physical parameters, a large amount of thermal phase mixing at T = Tc is 
expected to occur during the quark-hadron phase transition in the early universe, as well 
for the plasma produced in heavy- ion collisions [|T7|. For high enough temperatures and 



low enough cooling, large-amplitude thermal fluctuations will populate the new minimum 
at = 0ft in the range Tc < T < Ti. Although these fluctuations which are in the form 
of subcritical hadron bubbles will shrink and flnally disappear, there will always be some 
non-zero number density of hadron bubbles at a given temperature T. In this work, we 
study the equilibrium density distribution of subcritical hadron bubbles for a wide spectrum 
of very weak to very strong flrst order QCD phase transition, using the formalism devel- 
oped in Refs. [p!^-pO|. It is found that the density of subcritical hadron bubbles builds up 
faster as the transition becomes weaker, leading, in some cases, to complete phase mixing 
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at T = Tc- Further, using reasonable values for the surface tension and correlation length 
as obtained from lattice QCD calculations, we find that (although large) the amount of 
phase mixing remains below the percolation threshold. Therefore, the quark-hadron phase 
transition will begin with the nucleation of critical-size hadron bubbles from a supercooled 
and inhomogeneous background of quark-gluon plasma. Since the background contains sub- 
critical hadron bubbles, the homogeneous theory of nucleation needs to be modified. In Ref. 

, an approximate method was suggested to incorporate this inhomogeneity by modelling 
subcritical bubbles as Gaussian fluctuations, resulting in a large reduction in the nucleation 
barrier. Here, we will study inhomogeneous nucleation in the framework of homogeneous 
theory, but with a reduced nucleation barrier that accounts for the inhomogeneity of the 
medium. Finally, we also briefly discuss possible implications of inhomogeneous nucleation 
to relativistic heavy-ion collisions and cosmology. 

The paper is organized as follows. In the next section, we begin with the discussion of 
a quartic double-well potential used to describe the dynamics of a first-order quark-hadron 
phase transition. The parameters of the potential are obtained in terms of relevant physical 
quantities such as critical temperature, surface tension and correlation length. In section III, 
we estimate the equilibrium fraction of subcritical hadron bubbles from very weak to strong 
first-order phase transitions. We also estimate the reduction in the nucleation barrier by 
incorporating the presence of subcritical bubbles in the medium. Using this reduced barrier, 
we study nucleation and supercoohng in section IV. Finally, we present our conclusions in 
section V. 



II. PARAMETERIZATION OF THE EFFECTIVE POTENTIAL 

We consider a general form of the potential (or equivalently, the homogeneous part of 
the Helmholtz free energy density) to study the quark-hadron phase transition in terms of 
a real scalar order parameter given by 

V{^,T) = a{T)4>^ -bT(f)'' + c(P\ (1) 



where b and c are positive constants. Ignatius et. al. |21[] use this parameterization to 



describe the phase transition from a QGP (symmetric phase) to a hadron phase (broken 
symmetry). The meaning of is obvious for a symmetry-breaking transition, but the same 
description can be used if no symmetry is involved. The order parameter could then be 
related to energy or entropy density. The parameters a, b and c are determined in terms of 
surface tension (a), correlation length [^) and critical temperature (T^). The potential has 
two minima, one at 0g = and the other at cph = {'ibT + \/9b'^T'^ — 32ac) j 8c, which in our 
case will represent quark and hadron phases respectively. These phases are separated by a 
maximum defined by 0^ = (36T - VQ&^T^ - 32ac) /8c. At T = Tc, 

V(0„Tc) = n0,,rc) = O, (2) 

having the required degeneracy. The above condition yields, 

aiTc) = b^T^/4c, MTc) = bTc/2c and 0^(Tc) = bTc/Ac. (3) 

Using these relations, the barrier height at Tc can be written as 
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(4) 



Therefore, if the parameter c is kept fixed, b can be varied to characterize a wide spectrum 
of very weak to very strong first-order phase transitions. The transition is strong enough 
for large Vb and very weak or close to second order as H — * 0. In the following, we relate 
the parameters b and c to the surface tension and the correlation length in the quark phase. 
The surface tension can be defined as the one dimensional action given by. 



cr 



dx 




+ ^(0) 



(5) 



Under the thin-wall limit, A = |V^(0) — V{(j)h)\ 0, the surface tension can be expressed as 



22 



a 



V2b^T^ 
48 



(6) 



Similarly, the correlation length around the quark phase is obtained using = 1/ ^JV"{(j))\^=Q 
= l/J2a{T). At the critical temperature, using Eq. (|]), we get 



2c 



c) 



bTc 



c 



From Eqs. (||) and we get 



(7) 



(8) 



in terms of the values of a and ^„ at Tc- The barrier height VJ, can now be written as 



(9) 



Thus, the barrier height is proportional to the ratio cr/^q. The transition becomes very weak 
as a decreases and increases. Here, we fix = 0.5 fm at T = Tc and vary a to investigate 
phase transitions with different strengths. The temperature dependence of a is deduced by 
equating the depth of the second minimum with the the pressure difference AP between the 
two phases at all temperatures. This yields an equation 



AP = pf^-pq 

= ViO) - Vi(l)H) 



a{T) - bT(j)h + cf^ 



h 'Th 



(10) 



which is solved to get the parameter a(T), giving the temperature dependence of ^g. The 
surface tension will also have small temperature dependence which we ignore, as we are not 
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going too far from the critical temperature. Thus, we have parameterized the free-energy 
density in terms of the surface tension, correlation length, critical temperature and equation 
of state, which can be obtained from lattice QCD calculations. The bag equation of state 
which is a good depiction of the lattice results is used to calculate the quark/hadron pressure 
Pq/h as follows 

Pg = agT'^-B, ph = ahT^, (11) 

where B = {ag — ah) is the bag constant. The quark phase is assumed to consist of a 
massless gas of u and d quarks and gluons, while the hadron phase contains massless pions. 
Thus, the coefficients Qg and at are given by Qg = 3777^/90 and at = 377^/90. The critical 
temperature is taken as Tq = 160 MeV. 

Fig. 1 shows the plot of V{(f)) as a function of (p at three different temperatures for 
a typical value of a = 30 MeV/fm^ and ^g(Tc) = 0.5 fm. At T = Tc, the potential is 
degenerate with a large barrier that separates the two phases. Below Tc, the phase (f) = (ph 
has lower free-energy density, and the QGP phase becomes metastable. Above Tc, the 
potential has a metastable minima at (p = (ph (hadron phase) as long as T remains below Ti. 
The temperature Ti [at which (ph = (pm and 9b^T^ = 32a(Ti) c] can be obtained analytically 
by solving Eq. (|T0|) as. 



B 



1/4 

Tc. (12) 



It may be mentioned here that the dynamics of the phase transition has also been studied 
in Ref. [|17| using a different form of the potential which has been parameterized as a fourth 



order polynomial in the energy density [g. This form is unsuitable over a wide range of 
temperatures due to the persistence of metastability at much above and below Tc- 



III. MODEL FOR LARGE- AMPLITUDE FLUCTUATIONS 



We closely follow the work of Refs. [16,|19,20| to estimate the equilibrium density dis- 
tribution of subcritical hadron bubbles by modeling them as Gaussian fluctuations with 
amplitude (pA and radius R 



5Ae 



and 



(1 



— e 



(13) 



The amplitude (pA is the value of the field at the bubble's core away from the quark phase. 
For smooth interpolation between the two phases in the system, cpA > (pm- The free energy 
of a given configuration can then be found by using the general formula p2 |. 



-XV<P{r)f + V{<P{r)) 



Using Eq. (|T3p and Eq. (0) in Eq. (0) we get 

Fg_,h = Fh = ahR + PhR^ and Fh^g 
where ah, Ph, C(g and Pg are given by 



Fg = agR + /3g R^^ 



(14) 



(15) 
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ah 



an 



3v^ 



'3/2-2 n 



(16) 



and 



+ 



72 
4 



3/2, 



'1 ^ 3^2 



4v^ 
~9~ 



4 cvr 



3/2, 



(17) 



It may be mentioned here that ah{= oiq) is positive and is much greater than (3h{q)- 
Therefore, the free energy grows hnearly for small values of R. Further, hadron bubbles of 
all configurations will be subcritical as long as j3h{q) is positive. At T = Tc, both [3h and Pq 
are positive for all amplitudes. However, below Tc, I3h may become negative for some values 
of (pA- For such configurations, the free energy has a maximum at Rm = \f(yh/SPh and these 
bubbles are not strictly subcritical. The same is true for Pg above Tc- We thus restrict the 
amplitudes (j)A to the range where I3h(q) is positive. If not exactly the same, the limits of 
integration (pmin and (pmax for (pA are found to be quite close to 0^ and 0/j respectively. 



A. Equilibrium fraction of subcritical bubbles 

There will be fluctuations from quark to hadron phase and back. To obtain the number 
density ua of subcritical bubbles, we define the distribution function / = O'^ha/ dRd(j)A where 
/(i?, 0^, t)dRd(f)A is the number density of bubbles with radius between R and R + dR and 



amplitude between (pA and 0yi + dcpA at time t. It satisfies the Boltzmann equation |T6|,pO 



The first term on the RHS is the shrinking term. Here, |f | is the shrinking velocity, which 
we assume to be given by the velocity of sound (=1/ V^) in a massless gas. The second term 
is the nucleation term where G is the Gibbs distribution function defined as F = / dRd(j)G. 
Here Th is the nucleation rate per unit volume of subcritical bubbles from the quark phase 
to the hadron phase. Similarly Fg is the corresponding rate from the hadron phase to the 
quark phase. The factor 7 is defined as the fraction of volume in the hadron phase and is 
obtained by summing over subcritical bubbles of all amplitudes and radii within this phase. 
The Gibbs distribution function is defined as |T^, 



G^/q = AT^e-^'^/«(^''^^)/^, (19) 

where A is of O ~ 1 

If the equilibration time scale is smaller than the expansion time scale of the system, 
we can obtain the equilibrium number density of subcritical bubbles by solving Eq. (|18D 
with df /dt = 0. Since the early universe expands at a much slower rate pO| , pJ| , the above 
assumption is quite reasonable in the context of the cosmological QCD phase transition. 
However, QGP produced during heavy ion collision may expand at a faster rate as compared 
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to the early universe. In this case, it is possible that the density distribution of the subcritical 
bubbles will not attain full equilibrium. For simplicity, we will assume an equilibrium 
situation so that the present results on the fraction of subcritical bubbles and phase mixing 
can be considered as an upper limit. Using the boundary condition f{R oo) = 0, we get 
the equilibrium distribution given by 

f{R, (Pa, T) = (1 - 7) Ws{R, (Pa, T) - 7 Wt{R, <Pa, T), (20) 

where 

Ws{R,(pA,T) = {A/\v\)T^ / e-("'^^'+^'^^")/^cii?', 

JR 

Wt{R, (pA, T) = {A/\v\)T^ r e-^^'^'+^'^^^/^rfi?'. (21) 
The equilibrium fraction 7 of volume occupied by subcritical bubbles is given by, 

nax r ^max ^TT 



r<Pmax r^max z^TT 

Rm..)= / —R'f{R,(PA,T)dRd(PA, (22) 

J Rm 

which is solved to get 



where 



rRm ax ZJ-TT f, 

IsiT)= / —R'WsiT){R,(pA,T)dRd(PA. (24) 

Here, (pmin and (pmax define the range within which both and Pg are positive. Rmin is 
the smallest radius of the subcritical bubbles taken as ^q, the correlation length of the fluc- 
tuations. The R integration should be carried out over all bubbles with radii from Rmin = 
to Rmax = 00. For very weak transitions, both a and (3 are very small and the R integration 
may not have good convergence. However, we found that the 7 value is maximized when 
Rmax is about 3 to 4 fm. Therefore, we use i?maa;=3.5 fm. This is a reasonable choice as 
bubbles with i? ~ will be statistically dominant and larger fluctuations have larger free 
energy and are exponentially suppressed. 

Fig. 2 shows the plot of the subcritical hadron fraction 7 as a function of a at T = 
and at a fixed value of C,q(Tc) = 0.5 fm. The fraction 7 has been estimated (dashed curve) 
assuming that, for a degenerate potential, Gh — Gg, as in Ref. |jl6|. This assumption is valid 



only for the configuration for which (pA = (ph- However, when we include other configura- 
tions in the range (pmin to (pmax, the integral It turns out to be always higher than Is at Tq. 
Therefore, 7 obtained using Gh 7^ Gg is always lower than when the approximation Gh = Gg 
is used. In both cases, the value of 7 increases with decreasing a i.e. when the transition 
becomes weak. It may be mentioned here that as per lattice QCD calculations without dy- 
namical quarks 0, a lies between 2 MeV/fm^ and 10 MeV/fm^. There would be 15% to 30% 
phase mixing corresponding to these a values, which is still below the percolation threshold 
(7 < 0.3). If 7 > 0.3, the two phases will mix completely, the mean-field approximation 
for the potential breaks down, and the phase transition may proceed through percolation 
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p!5| , |20[] . However, for a surface tension in the range 2 MeV/fm^ < < 10 MeV/fm^, the 
phase transition will proceed through the formation of critical-size hadron bubbles from a 
supercooled metastable QGP phase. Since the QGP phase is no longer homogeneous, the 
dynamics of the phase transition will be quite different from what is expected on the basis 
of homogeneous nucleation theory ||T6|. We refer to it as "inhomogeneous nucleation." 

We would also like to mention here that the present results are in disagreement with 
the findings of Ref. []T^, where a large fraction of subcritical hadron phase was found at 
and above Tq- This scenario is highly unrealistic and probably could be due to the choice 
of the potential parameterization, which shows a metastable hadron phase much above Tc- 
Therefore, the authors of Ref. [|1^] found a finite fraction of hadron phase at temperatures as 
high as twice Tc- Furthermore, the value of 7 strongly depends on how the shrinking term 
is incorporated in the calculation. In our case, it is proportional to the gradient {Of /OR) 
that appears in the kinetic equation (|18D in a natural way, whereas in Ref. [1^, a specific 
assumption is made to take into account the shrinking of the hadronic volume. 



B. The total free energy of subcritical bubbles and the nucleation barrier 

The nucleation rate in the standard theory which neglects phase mixing, is given by 

/~Ar^e-^c/T_ (25) 

Here Fq is free energy needed to form a critical bubble in the homogeneous metastable back- 
ground. For an arbitrary thin-wall spherical bubble of radius R and amplitude (pthin ^ 4'hi 
the free energy of the bubble takes the well-known form 

477- 

Fthin(i?) = -—R:'/W + 47ri?V. (26) 

In the above, AV is defined as the difference in free-energy density between the background 
medium and the bubble's interior. For a homogeneous background (metastable), we can 
write, 

AV = AVo = V{0) - VM. (27) 

If there is significant phase mixing in the background metastable state, its free energy is 
no longer ^(0). One must also account for the free energy density of the nonperturbative 



large amplitude fluctuations. Following Ref. |16|, we write the free energy density of the 



metastable state as 1^(0) + J-'sc, where J-'sc is the extra free energy density which can be 
estimated from the density distribution of subcritical bubbles as follows: 



^sc~/, / Fh{R,(j)A,T) f{R,(j)A,T)dRd^A, 

rmn 111 

max r F^max C'pmax f F^max 



7) / / Ff,WsdRd(j)A-l / / Fj,WTdRd(j)A. (28) 



R„ 



Once we know the hadronic fraction 7 and the free energy for a bubble of a given radius 
R and amplitude (pA, we can estimate the free-energy density correction due to the presence 
of Gaussian subcritical bubbles. 
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Since, for a critical size bubble, dF/dR\R^ = 0, we can use Eq. (p6D to obtain the free 
energy needed to form a thin-wall critical bubble in a background of subcritical bubbles, 

47r , 2(T , , 

Fc = —aRl, Rc = -TTT (29) 

For a very strong first-order phase transition, the subcritical bubbles are suppressed {J^^c 
0), and both Fq and Rc approach the homogeneous background expression. However, in the 
presence of subcritical bubbles, extra free energy becomes available in the medium, reducing 
the nucleation barrier. In other words, the extra background energy enhances the nucleation 
of critical bubbles. To illustrate this, we have plotted Fq/T and 7 as a function of T /Tq in 
Figs. 3 to 5 with a values of 50 MeV/fm^, 30 MeV/fm^ and 10 MeV/fm^, respectively, which 
are widely used in the literature. As evident, with decreasing temperature, the nucleation 
barrier decreases and the subcritical hadron fraction 7 increases. The reduction in barrier 
height due to J^sc (or due to 7 ) is more significant for lower values of a, corresponding to a 
weaker transition. Since the height of the nucletaion barrier decreases, the nucleation rate 
will also be enhanced, reducing the amount of supercooling further. The time evolution of 
the temperature and the supercooling are discussed in the next section. 



IV. NUCLEATION AND SUPERCOOLING 

As mentioned before, the background metastable state is inhomogeneous due to subcrit- 
ical hadron bubbles. It is now possible to study the kinetics of the nucleation of critical 
hadron bubbles using the corrected nucleation rate, as obtained in the previous section. In 
the present work, the prefactor in the nucleation rate is taken as AT'^ [see Eq. (P3|)]. In 
our previous work, 0, we have used a prefactor derived by Csernai and Kapusta for 
a dissipative QGP. In Ref. |]23|, Ruggeri and Friedman had derived a prefactor for a non- 
dissipative QGP. Recently, using a more general formalism, we have also derived a prefactor 
^ which has both dissipative and non-dissipative components corresponding to Ref. |^ 



and Ref. respectively. However, for consistency with the subcritical formalism, we use 



a more generic form Jq = AT^, with A a constant of order unity, as used in many studies 
of quark-hadron phase transition (see, for example, Refs. [|I^,|r^). The question of how to 
estimate the prefactor appearing in the nucleation rate of subcritical bubbles remais open. 
Using the nucleation rate /(T), the fraction h of space which has been converted to hadron 
phase due to nucleation of critical bubbles and their growth can be calculated. If the system 
cools to Tc at a proper time Tc, then at some later time r the fraction h is given by 

h{T) = r dr'I (T(r')) [1 - h{r')] V{t', t). (30) 

Jtc 

Here, V{t' , r) is the volume of a critical bubble at time r which was nucleated at an earlier 
time r'; this takes into account the bubble growth. The factor [1 — h^r)] accounts for the 
available space for new bubbles to nucleate. The model for bubble growth is simply taken 



as 25 



V{t',t) = ^ (^RciTir')) + f^^dr"v{T{r")))\ (31) 
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where v{T) = 3[1 — T/Tc]^/^ is the velocity of the bubble growth at temperature T P,|26 
The evolution of the energy density in 1+1 dimensions is given by 

de oj , , 

d7 + 7 = ''- (^2' 

The energy density e, enthalpy density u and the pressure p in pure QGP and hadron phases 
are given by the bag model equation of state. In the transition region, the e and at a time 
r can be written in terms of the hadronic fraction as 

e(r)=e,(T) + [e,(T)-e,(T)]Mr), 
uj{t) = oOgiT) + K(T) - LUg{T)]h{T). (33) 

Equations (|30|), (p2D, and ( ^3]) are solved to get the temperature as a function of time in the 
mixed phase with the initial conditions for temperature Tq = 250 MeV and proper time 
To = 1 fm/c at T(7=160 MeV. After getting T and /i as a function of time, the density of 
nucleating bubbles at a time r can be obtained in our model as 

N{t)= r dr'I{T{r'))[l-h{r')]. (34) 

Jtc 

The density would increase as the temperature drops below Tc and would ultimately 
saturate as h increases. 

Figure 6 shows the temperature variation as a function of proper time r at a = 50 
MeV/fm^. As the system cools below T^, the nucleation barrier decreases and also 7 in- 
creases. If only homogeneous nucleation (dashed curve) is considered, the system will su- 
percool up to 0.945 Tc- At this temperature, the hadronic fraction 7 has reached 10 % (See 
Fig. 3), which corrects the amount of supercooling (solid curve) by about ~10 % (up to 
0.95 Tc). Figure 7 shows a similar study at cr = 30 MeV/fm^. Since the nucleation barrier 
reduces with decreasing cr, the system supercools only up to 0.98 Tc under homogeneous 
nucleation. The hadronic fraction 7 corresponding to this value is ~ 12 — 13 % (See Fig. 4) 
which reduces the amount of supercooling by about ~ 20% (up to 0.984 Tc). For a around 
10 MeV/fm^, the supercooling will be reduced further (up to ~ 0.997Tc). Lattice QCD 
calculations predict a surface tension even smaller than 10 MeV/fm^, indicating a very weak 
first order transition. Although supercooling will be reduced further with decreasing a, we 
do not use very small a due to increased numerical inaccuracy. Further, it may be mentioned 
here that, although the fraction 7 grows with decreasing a, we never encountered 7 greater 
than 0.3: we remained within the sub-percolation regime throughout our analysis. 

Apart from a, the amount of supercooling also depends on Tc, the time taken by the 
system to cool from Tq to Tc- In QGP phase, the solution of Eq. (|32|) (T^r =constant) 
predicts Tc = tq{TqITcY . The choice of tq = Ifm, To=250 MeV and Tc=160 MeV results 
in rc=3.8 fm/c. However, formation of QGP with higher initial temperature (as high as 
3 to 4 times Tc resulting in large r^) can not be ruled out at RHIC and LHC energies 
p7| . Therefore, we have also studied the effect of Tc on supercooling, specifically, on the 
hadronization rate as well as on the density of nucleating bubbles. Figs. 8(a) and 8(b) show 
the plots of N{t) and h{T) as a function of r for two typical values of (3.8 fm/c and 25 
fm/c) corresponding to cr=10 MeV//m^ both with (solid curve) and without (dashed curve) 
inhomogeneity corrections. A general observation (both with and without correction) is that 
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the amount of supercooling, the rates of hadronization and bubble nucleation are reduced 
when Tc becomes larger. Although supercooling reduces with increasing Tc, the system will 
get reheated at an earlier temperature and also will encounter a larger nucleation barrier 
as compared to the case when is small. As a result, the rate of hadronization and also 
the rate of increase of density of the nucleating hadron bubbles will proceed at a slower 
rate when Tc is large. However, the reverse happens when the inhomogeneity correction is 
applied. Even though the medium gets heated up earlier, the reduction in the barrier height 
is quite significant as T approaches Tq. Another parameter that affects both N{t) and /i(r) 
is the expansion rate of the medium, i.e., the rate of change of temperature between and 
r, which also depends on r^. The overall effect is that both N{t) and h{T) rise faster as 
compared to their homogeneous counterparts (see Fig. 8 for Tc=25 fm/c), particularly when 
Tc is very large. (Compare the left and right curves on Fig. 8.) The increase in rates of 
N{t) and /i(r) is also larger for small a at large r^. 

For weak enough transition, the presence of inhomogeneity may also affect several observ- 
ables which can be detected experimentally. For example, the faster rate of hadronization 
at large Tc as compared to its homogeneous counterpart will lower the amount of entropy 
production, which, in turn, will affect the final hadron multiplicity distributions. Although 
not studied here, the bubble size distribution will also be affected by the dynamics of nu- 
cleation p. Since the nucleating bubble will act as a source of pion emission, the effect of 
inhomogeneity can also be inferred through interferometry measurements. 

In a cosmological context, the value of Tc is much larger than what was quoted here. 
Since the presence of inhomogeneities weakens the transition, more critical bubbles will be 
nucleated per unit volume, decreasing the inter-bubble distance, {d ~ A^~^/^); the pres- 
ence of subcritical bubbles can be thought as seeds for nucleation. As a consequence, the 
transition will produce smaller fluctuations in baryon number, protecting homogeneous nu- 
cleosynthesis. Although the present study is indicative enough of the reduction in the mean 
inter-bubble separation as compared to homogeneous nucleation, a quantitative estimate 
would require a more detailed analysis, including expansion. However, since the cosmologi- 
cal expansion rate is typically much slower than the subcritical bubble nucleation rate, we 
believe our results for the inter-bubble distance will carry on in this case as well. 

V. CONCLUSIONS 

We have estimated the amount of phase mixing due to subcritical hadron bubbles from 
very weak to very strong first-order phase transitions. With a reasonable set of values for 
the surface tension and correlation length (as obtained from lattice QCD calculations), we 
found that phase mixing is small at T = Tc , building up as the temperature drops further. 
We have shown that the system does not mix beyond the percolation threshold, allowing 
us to describe the dynamics of the phase transition on the basis of homogeneous nucleation 
theory with a reduced nucleation barrier. Accordingly, we have found an enhancement in 
the nucleation rate which further reduces the amount of supercooling. Although we have 
not included cosmological expansion in our analysis, we believe that our results indicate 
that the presence of an inhomogeneous background of subcritical bubbles will decrease the 
inter-bubble mean distance, and thus the fluctuations in baryon number which could damage 
homogeneous nucleosynthesis. 



11 



We have assumed that the equihbration time-scale for sub critical fluctuations is much 
larger than the cooling time-scale of the system. This may be the case for a quark-hadron 
phase transition in the early universe, where the expansion rate is quite slow. In the case of 
QGP produced at RHIC and LHC, the coohng rate is much faster than cosmological time- 
scales, and the subcritical bubbles density distribution may not attain full equilibration. We 
are presently investigating this issue in more detail. However, the present results should 
provide an upper bound on the fraction of subcritical hadron bubbles and their effect on 
supercooling and nucleation rates. 
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Figure Captions 

Fig. 1 The effective potential as a function of order parameter at, below and above T^. 

Fig. 2 Subcritical hadronic fraction 7 as a function of surface tension a. 

Fig. 3 The nucleation barrier Fc/T for critical bubbles with (solid line) and without 
subcritical bubble correction (dashed curve) as function of temperature for a — 50 MeV/fm^ 
is shown in upper panel. Corresponding subcritical hadron fraction 7 is shown in the lower 
panel. 

Fig. 4 Same as Fig. 3 but at a = 30 MeV/fm^. 
Fig. 5 Same as Fig. 3 but at tr = 10 MeV/fm^. 

Fig. 6 The temperature variation as a function of proper time with (solid curve) and 
without subcritical bubble correction (dashed curve) for a— 50 MeV/fm^. 
Fig. 7 Same as Fig. 6 but at a= 30 MeV/fm^. 

Fig. 8 (a) Density of nucleating bubbbles as a function of proper time with (solid 
curve) and without subcritical bubble correction (dashed curve) for a— 10 MeV/fm^ (b) 
the hadronic fraction 7 as a function of r. The left curves are for tc — 3.8 fm/c and the 
right curves for tc — 25 fm/c. 
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